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SUMMARY 


A chemical kinetic reaction mechanism for the combustion of silane/hydrogen was 
assembled through the analogy of silane to methane. The mechanism was refined by 
comparing the observed kinetic behavior as determined in shock tube studies with that 
predicted by the mechanism. Good agreement between the experimental and the calcu- 
lated results was achieved by the adjustment of the rate coefficient for the reaction 
SiH.j + C> 2 -* SiH 2 0 + OH. The reaction mechanism was used to theoretically investigate 
the ignition characteristics of si lane /hydrogen mixtures. The results revealed that 
over the entire range of temperature examined (800 K to 1200 K), substantial reduc- 
tion in ignition delay time is obtained when silane is added to hydrogen. 


INTRODUCTION 

Much attention is being paid to the use of silane (SiH^) as an ignition aid for 
hydrogen-fueled supersonic combustion ramjet engines. Experiments in the Langley 
Hypersonic Propulsion Test Cells and at contractor facilities have positively demon- 
strated that the addition of silane to the hydrogen fuel enhances combustion and 
helps sustain stable combustion (Beach et al. 1980). In addition, recent shock tube 
studies (McLain et al. 1983) have shown that adding silane to hydrogen-oxygen- 
nitrogen mixtures decreases the ignition delay time substantially. To take full 
advantage of the effect of silane, a better understanding of the chemical reaction 
process by which silane participates in the combustion process is needed. A chemical 
kinetic mechanism which describes the essential features of the silane/hydrogen com- 
bustion process would be a useful tool not only for assessing the sensitization capa- 
bilities at other conditions but also for the analysis and interpretation of combus- 
tor results . 

Beach et al. (1980) assembled a preliminary chemical kinetic mechanism for 
silane oxidation obtained by direct analogy with the methane oxidation reaction mech- 
anism. The mechanism was tested by comparing the calculated consumption of silane 
with experimental results from a low temperature (500 K to 700 K) isothermal flow 
reactor. Rate coefficients for five reactions were adjusted to obtain agreement 
between theoretical and experimental results. Since a large number of rate coeffi- 
cient adjustments and some unusually large values for some of the rate coefficients 
were required to obtain agreement, it is thought that heterogeneous processes as well 
as gas phase homogeneous processes were occurring in the flow reactor. In this 
report, a detailed chemical mechanism for the oxidation of si lane /hydrogen mixtures 
is presented and discussed. Rate coefficients in this mechanism are refined by com- 
paring calculated ignition delay times with the experimental data reported by McLain 
et al. (1983) . 


SILANE/HYDROGEN REACTION MECHANISM 

Since very little kinetic and mechanistic information is available on the high 
temperature silane oxidation process, the basic features were taken from the methane 
oxidation mechanism. This reaction has been studied for many years, and the essen- 
tial reaction paths, intermediates, and products are well known. The silane molecule 
has a structure similar to that of the methane molecule, and the silane oxidation 



products (SiO, SiO,,) are similar to the methane oxidation products (CO, C0 2 ) . There- 
fore, it seems reasonable to assume that similar reactions and reaction intermediates 
would occur. Because silicon bonding and carbon bonding are not identical, it is 
recognized that reactivity and the nature of the intermediate species may be differ- 
ent from the methane scheme. It was not the intent, however, to assemble an exten- 
sive list of reactions such as those recently developed for methane (Westbrook 1979; 
Olson and Gardiner 1978). Instead, the objective was to develop a silane mechanism 
which contained enough detail to describe the essential features of the oxidation 
process. The mechanisms used to describe the combustion of silane and hydrogen are 
shown in tables I and II. Rate coefficients for the hydrogen mechanism were selected 
from the literature. Whenever possible, the same was done for silane. However, 
since very little kinetic information on silane is available, most rate coefficients 
were estimated. 


TABLE I.- SILANE REACTION MECHANISM 


Reaction 


SiH 4 a SiH 2 + H 2 
SiH 4 + 0 2 > SiH 3 + H0 2 
SiH 4 + H0 2 a SiH 3 + H 2 0 2 
H + SiH 4 a h 2 + SiH 3 
0 + SiH 4 a OH + SiH 3 
OH + SiH 4 a h 2 0 + SiH 3 
H + SiH 3 a siH 2 + H 2 
0 + SiH 3 a- SiH 2 0 + H 
OH + SiH 3 a SiH 2 0 + H 2 
SiH 3 + 0 2 a SiH 2 0 + OH 
SiH 2 + 0 2 a- HSiO + OH 
H + SiH 2 o •> h 2 + HSiO 
0 + SiH 2 0 a- OH + HSiO 
OH 4* SiH 2 0 a- H 2 0 + HSiO 
H + HSiO a H 2 + SiO 
0 + HSiO a OH + SiO 
OH + HSiO a H 2 0 + SiO 
HSiO + + M a H + SiO + M 
HSiO + 0 2 a SiO + H0 2 
SiH 2 0 + H0 2 a HSiO + H 2 0 2 
SiO + 0 + M a Si0 2 + M 
SiO + OH a Si0 2 + H 
SiO + 0 2 a Si0 2 + 0 


Rate 


coefficient 


* 


6.0 

X 

10 13 

2.0 

X 

10 1 1 

3.0 

X 

10 1 2 

1.5 

X 

10 13 

4.2 

X 

10 12 

8.4 

X 

10 12 

1 .5 

X 

10 1 3 

1 .3 

X 

10 14 

5.0 

X 

10 1 2 

8.6 

X 

10 14 

1 .0 

X 

10 14 

3.3 

X 

10 14 

1 .8 

X 

10 1 3 

7.5 

X 

1 o 1 2 

2.0 

X 

10 1 4 

1 .0 

X 

10 1 4 

1 .0 

X 

10 1 4 

5.0 

X 

10 1 4 

3.0 

X 

10 1 2 

1 .0 

X 

10 1 2 

2.5 

X 

10 15 

4.0 

X 

10 1 2 

1 .0 

X 

10 1 3 


exp( -54960/RT) 
exp( -44000/RT) 
exp( -5600/RT) 
exp( -2500/RT) 
exp(-1 600/RT) 
exp(-1 OO/RT) 
exp ( -2500/RT) 
exp( -2000/RT) 

exp(-1 1400/RT) 
exp( -3700/RT) 
exp( -1 0500/RT) 
exp( -3080/RT) 
exp ( -1 70/RT) 


exp ( -29000/RT) 

exp( -8000/RT) 
exp( -4370/RT) 
exp( -5700/RT) 
exp( -6500/RT) 


The rate coefficient is defined by k = A exp(-E/RT) 
where A is in s - ”' for unimolecular reactions, cm 3 /mole-s 
for bimolecular reactions, and cm /mole -s for termolecular 
reactions. The activation energy E is in cal/mole, 

R = 1 .987 cal/mole-K, and T is temperature in K. 

^Determined in this study. 

*M is a third-body molecule. 
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TABLE II.- HYDROGEN REACTION MECHANISM 


Reaction 

Rate coefficient* 

Reference 

H 2 + ° 2 -► OH + OH 

1.7 x 10 13 exp (-481 OO/RT) 

Jachimowski and Houghton 1971 

H + 0 2 * OH + O 

1.2 x 10 17 T -0 * 91 exp( -1 6600/RT) 

Schott 1973 

O + H 2 ■* OH + H 

2.1 x 10 14 exp(-1 3750/RT) 

Schott 1973 

OH + H 2 -> H 2 0 + H 

3.2 x 10 7 T 1,8 exp( -3030/RT) 

Gardiner et al. 1974 

OH + OH -)■ H 2 0 + O 

5.5 x 10 13 exp( -7000/RT) 

Gardiner et al. 1973 

H + OH + + M > H 2 0 + M 



M = HO 

1.4x1 0 23 T -2,0 

Baulch et al. 1972 

M = Ail other 

2.2 x 10 22 T -2,0 

Baulch et al. 1972 

H+H + M->-H 9 + M 



M = h 2 

2.6 x 1 0 1 2 T -1 *° 

Baulch et al. 1972 

M = H 2 0 

3.9 x 10 12 T" 1 *° 

Baulch et al. 1972 

M = All other 

6.5 x 10 17 T -1 *° 

Baulch et al. 1972 

H + 0 2 + M + H0 2 + M 



M = h 2 

7.0 x 1 o 1 8 T -1 *° 

Slack 1977 

M = H 2 0 

4.4 x 10 19 T" 1 *° 

Slack 1977 

M = All other 

3.2 x 10 18 T -1 *° 

Slack 1977 

OH + H0 2 > H 2 0 + 0 2 

5.0 x io 13 exp( -1 000/RT) 

Lloyd 1974 

H + H0 2 ->■ H 2 + 0 2 

2.5 x io 13 exp ( -700/RT) 

Lloyd 1974 

H + H0 2 ->• OH + OH 

2.0 x io 14 exp ( -1 900/RT) 

Lloyd 1974 

0 + H0 2 -> OH + 0 2 

5.0 x io 13 exp(-1 000/RT) 

Lloyd 1974 

ho 2 + ho 2 -*■ h 2 o 2 + o 2 

2.0 x io 12 

Slack 1977 

H0 2 + H 2 ■> H 2 0 2 + H 

3.0 x 10 11 exp( -1 8 700/RT) 

Slack 1977 

OH + H 2 0 2 ■> H0 2 + H 2 0 

1.0 x io 13 exp( -1 800/RT) 

Baulch et al. 1972 

M + H 2 0 2 -> OH + OH + M 

1.2 x 10 17 exp( -45500/RT) 

Baulch et al. 1972 


The rate coefficient is defined by k = AT 11 exp(-E/RT). The rate coefficient units 
are s - ^ for unimolecular reactions, cm 3 /mole-s for bimolecular reactions, and cm®/mole^-s 
for termolecular reactions. Hie activation energy E is in cal/mole, R = 1.987 cal/mole-K, 
T is temperature in K, and n is a temperature exponent. 

*M is a third-body molecule. 
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When gaseous mixtures of reactants are suddenly heated by a shock wave, ignition 
is preceded by a period during which the concentration of free radicals grows very 
rapidly. The initial phase of this ignition period is often characterized by decom- 
position of the fuel molecules. Whereas the high temperature decomposition of 
methane occurs primarily through the rupture of a C-H bond to form a methyl radical 
(CH 3 ) and a hydrogen atom, the initial reaction in the decomposition of silane is the 
molecular elimination of hydrogen 


SiH 4 ->■ SiH 2 + H 2 

Newman et al. (1979) investigated the homogeneous gas phase decomposition of 
silane using the single-pulse shock tube technique. They found decomposition rates 
which were in good agreement with the high pressure unimolecular rate coefficient, 
k = 3 x 10 exp ( -59600/RT) , s - ^ . In table I, the decomposition rate coefficient is 

CO 

given for the temperatures and pressures of the si lane /hydrogen shock tube studies 
reported by McLain et al. (1983). This rate coefficient was obtained from k and 
the results of the RRKM unimolecular reaction theory calculations given by Newman 
et al. (1979) . 

Also included in the silane mechanism as an initiation step is the hydrogen 
abstraction reaction 


SiH 4 + 0 2 a- SiH 3 + H0 2 


The rate coefficient for this reaction was estimated by setting the activation energy 
equal to the heat of reaction AH at 298 K and using the preexponential factor 
reported for the analogous CH 4 reaction (Tsuboi 1976). The rate coefficient for the 
reaction 


SiH 4 + H0 2 ■> SiH 3 + H 2 0 2 


was estimated by using the relation E = AH + E(reverse) to calculate the activation 
energy; the preexponential factor was set equal to the preexponential reported for 
the analogous CH 4 reaction (Westbrook 1979). In the above relation, AH is the heat 
of reaction at 298 K and E(reverse) is the activation energy for the reverse 
(exothermic) reaction. The latter was calculated using Hirschfelder ' s rule (Wayne 
1969); E(reverse) = 0.055 Dq_ h , where D Q _ H is the OH bond energy in the H 2 0 2 mole- 
cule. The rate coefficients for the primary chain propagating and branching reac- 
tions in the silane mechanism 


H + SiH 4 -► H 2 + SiH 3 
O + SiH 4 -> OH + SiH 3 
OH + SiH 4 -> H 2 0 + SiH 3 
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were obtained from the studies reported by Atkinson and Pitts (1978) and Arthur and 
Bell (1978) for temperatures between 297 K and 438 K. The rate coefficient for the 
reaction 


H + SiH 3 ->• H 2 + SiH 2 


was set equal to the rate coefficient for the H + SiH^ reaction. Olson and Gardiner 
(1978) made the analogous assumption when developing a detailed methane mechanism. 

The rate coefficients for the SiH^-consuming reactions 


SiH 3 + O ■> SiHjO + H 
SiH 3 + OH -> SiH 2 0 + H 2 
SiH 3 + 0 2 -»• SiH 2 0 + OH 


were set to the rate coefficients for the analogous CH 3 reactions as reported by 
Peeters and Mahnen (1973a), Peeters and Vinckier (1974), and Jachimowski (1974), 
respectively. The rate coefficient for the SiH 2 -consuming reaction 


SiH 2 + 0 2 -»■ HSiO + OH 


was set equal to the rate coefficient for the analogous CH 2 reaction as reported by 
Peeters and Mahnen (1973b). 

The rate coefficients for the reactions 


H + SiH 2 0 -+ H 2 + HSiO 
0 + SiH 2 0 ■+ OH + HSiO 
OH + SiH 2 0 ■> H 2 0 + HSiO 
H + HSiO -> H 2 + SiO 
0 + HSiO -> OH + SiO 
OH + HSiO -> H 2 0 + SiO 
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were set equal to the rate coefficients for the analogous Cf^O/HCO reactions as 
reported by Dean et al. (1980) in their study of the oxidation of CE^O behind shock 
waves. The rate coefficient for the decomposition reaction 


HSiO + M -> SiO + H + M 


was estimated by setting the activation energy equal to the heat of reaction AH at 
298 K and the preexponential factor for the analogous HCO reaction (Jachimowski 
1977) . 

The rate coefficients for the reactions 


HSiO 0 2 ^ SiO 4- HO 2 
SiH 2 Q + H0 2 > HSiO + H 2 0 2 


SiO + O + M + Si0 2 + M 


were set equal to the rate coefficients for the analogous carbon reactions by using 
the expressions recommended by Olson and Gardiner (1978), Jensen and Jones (1978), 
and Lloyd (1974), respectively. The rate coefficients for the Si0 2 -f orming reactions 


SiO + OH -> Si0 2 + H 

sio + o 2 -* sio 2 + 0 


were estimated using Hirschfelder 1 s rule to calculate the activation energies. The 
preexponential factors were set equal to the preexponential factors for the analogous 
CO reactions as reported in Gardiner et al. (1973) and Jensen and Jones (1978). It 
should be noted that SiO is not a stable species like CO, and that its reactivity is 
likely to be different. 

The hydrogen reaction mechanism listed in table II was taken from Slack (1977). 
This mechanism successfully reproduced experimental ignition time data from shock 
tube studies of hydrogen-oxygen mixtures over the temperature range 980 K to 1176 K. 


MECHANISM EVALUATION AND REFINEMENT 

As previously noted, the experimental shock tube results reported by McLain et 
al. (1983) were used to evaluate and refine the silane mechanism. McLain et al. 
(1983) reported ignition delay times for two silane/hydrogen mixtures, 2 percent 
SiH 4 , 8 percent H 2 , 4 percent 0 2 , and 86 percent N 2 (Overall equivalence ratio = 2), 
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and 1.68 percent SiH^, 6.72 percent H 2 , 6.74 percent 0 2 , and 84.86 percent N 2 
(Overall equivalence ratio = 1), at 1.25 atm and 1.35 atm, respectively, for the 
temperature range 800 K to 1050 K. Numerical simulation of the shock tube experi- 
ments was carried out using the chemical kinetic computer program described by McLain 
and Rao (1976). The program was operated in a constant volume mode to simulate the 
conditions behind a reflected shock wave. The ignition delay time is defined as it 
is in the shock tube experiments; that is, it is defined to be the elapsed time 
between the heating of the gas mixture by the reflected shock wave and the sudden 
pressure increase due to combustion. The rate coefficients for the reverse reactions 
were calculated within the program by using the forward rate coefficient and the 
appropriate thermochemical data. 

Thermochemical data for the hydrogen, oxygen, and nitrogen species were taken 
from the JANAF tables (1971) and were used in the analytical form described by 
Wakelyn and McLain (1975). Thermochemical data for the species SiH siH, SiO, and 
Si0 2 were obtained from the tables prepared by McBride et al. (1963)* The thermo- 
chemical data for the other silicon-containing species which appear in the silane 
mechanism were not available and had to be estimated. 


The thermochemical data for the radicals SiH^ and SiH 2 were obtained by assuming 
that their heat capacities at constant pressure were equal to the heat capacities for 
the analogous carbon-containing species, CH 3 and CH 2 ( Bahn 1973). The enthalpies of 
formation at 298 K were determined usinq the bond dissociation energies, D H ^g^_ H an< ^ 

Djj reported by Doncaster and Walsh (1981). The standard entropies at 

298 K (S°) were determined using the relation 


S°(SiH x ) = S°(SiH x _ 1 ) + S°(H) - AS° 


where AS° is the entropy change for the chemical change CH^ -> CH X _ 1 + H. In other 
words, the entropy changes for the reactions 


SiH 3 -> SiH 2 + H 
SiH 2 -y SiH + H 


were assumed to be equal to the entropy changes for the analogous carbon reactions 


CH 3 ->• CH 2 + H 
CH 2 -y CH + H 
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Thermochemical data for Sil^O and HSio were obtained by assuming that the heat 
capacities are the same as for the analogous carbon -containing species CH 2 0 and HCO. 
The enthalpies of formation at 298 K, AH°, were estimated using the relations 


AHj(SiH 2 0) = AHf(HSiO) + AHf(H) - %-HSiO 


AH f ( HSiO) = AH f (SiO) + AH f (H) - D H _ si0 


where the bond dissociation energies D H _ HSi0 and D H _ gi0 were set equal to D H _ HC0 and 
D u respectively. The standard entropies at 298 K were estimated using the 

ri"“UU , 

relations 


S°(SiH 2 0) = S° ( HSiO) + S°(H) - AS' 

S° ( HSiO) = S° ( SiO) + S°(H) - AS" 

where AS' and AS" are the standard entropy changes for the reactions 
CH 2 0 > HCO + H 
HCO + H + CO 


All the estimated thermochemical data were put into the polynomial format required by 
the chemical kinetic computer program. 

The results of the computer simulations of the shock tube experiments for the 
two si lane /hydrogen mixtures are plotted in figures 1(a) and 1(b), together with the 
corresponding experimental results. The results which were calculated with the esti- 
mated rate coefficients show ignition delay times much longer than the observed 
times. As a first step toward refining the model, a sensitivity study was carried 
out to determine which reactions had the largest influence on the calculated ignition 
times. The sensitivity study was performed by individually varying the rate coef- 
ficients for all the reactions in the silane mechanism so that they had values 
between 0.1 and 10.0 times the values listed in table I. 

The results of the sensitivity study were virtually identical to similar sensi- 
tivity studies carried out with the methane oxidation mechanism (Jachimowski 1974; 
Bowman 1975). This could have been expected since the silane reaction mechanism was 
patterned after the methane oxidation scheme, and many of the rate coefficients are 
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Ignition delay time, us Ignition delay time, us 



2 percent SiH^, 8 percent H 2 , 4 percent 0 2 , and 
86 percent N 2 mixture; Pressure = 1 .25 atm. 



(b) 1.68 percent SiH^, 6.72 percent H 2 , 6.74 percent O 2 , 
84.86 percent mixture; Pressure = 1.35 atm. 


Figure 1 .- Comparison between calculated and measured ignition 

delay times. 



similar. Varying the rates of the following reactions had the largest influence on 
the calculated ignition delay times. 


SiH 4 a- SiH 2 + H 2 
SiHg + 0 2 -> SiH 2 0 + OH 


SiO + OH •* Si0 2 + H 


The results of the sensitivity study did not indicate how much each of these 
rate coefficients should be modified to bring calculated results into agreement with 
observed ignition delay times. Good agreement could be obtained by simultaneously 
increasing the rate coefficients for all three reactions or by adjusting the rate 
coefficient for any one of the three reactions. However, some of the adjustments did 
not yield reasonable rate coefficient expressions. For example, the adjustment 
required to the rate of the decomposition reaction 


SiH 4 SiH 2 + H 2 


resulted in a rate coefficient several times the high pressure (upper limit) expres- 
sion reported by Newman et al. (1979). Similarly, the adjustment required to the 
rate of the reaction 


OH + SiO a- Si0 2 + H 


resulted in a preexponential factor greater than 5 x 10 15 , which is an unreasonably 
large value for a bimolecular reaction. Consequently, it was decided that it would 
be more prudent to modify the rate coefficient for the reaction 


SiH 3 + 0 2 SiH 2 0 + OH 


which had the largest influence on the calculated ignition delay times. The rate 
coefficient which gave the best agreement between the calculated and experimental 
results was k = 8.6 x 10*^ exp( -1 1 400/RT) , cm^/mole-s. The calculated results for 
the two mixtures using the modified rate coefficient are plotted in figures 1(a) 
and 1 (b) . The agreement between the calculated and observed ignition delay times is 
very good. 

Because of the scarcity of chemical kinetic data at high temperatures, it cannot 
be claimed that the mechanism presented here is unique. It is possible that other 
reactions that were not included in the mechanism may be important. Still, the mech- 
anism does accurately reproduce experimental results for silane/hydrogen oxidation in 
a shock tube with only a minimal amount of reaction rate modification; this suggests 
that the proposed mechanism and assigned rate coefficients are reasonably correct. 
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Nonetheless, experimental verification of the individual reactions and rate coef- 
ficients is still required. 


IGNITION CHARACTERISTICS OF SILANE/ HYDROGEN MIXTURES 

To acquire a better understanding of the effect of silane on the ignition char- 
acteristics of hydrogen-air mixtures, a series of calculations were made for a range 
of si lane /hydrogen mixtures, from pure hydrogen to pure silane. For each mixture, 
ignition delay period calculations were performed at several initial temperatures 
over the range 800 K to 1200 K for a constant pressure of 1 atm. The ignition delay 
time was defined as the time to achieve ignition as signified by the sudden increase 
in temperature. This definition gave ignition delay times which are essentially the 
same as those defined to be the time required to reach a temperature equal to the 
initial temperature plus 5 percent of the total temperature rise from the initial 
state to the final state (Beach et al. 1980). The overall fuel-air equivalence 
ratios ( 4>) considered were 1.0 and 0.5. The si lane /hydrogen mechanism with the modi- 
fied rate coefficient was used. 

The results of these calculations are given in figures 2(a) and 2(b) . The 
sensitization effect of silane is shown to be somewhat greater in stoichiometric 



(a) <|> = 1 . (b) <|> = 0.5. 

Figure 2.- Ignition delay time as a function of temperature for various 

si lane /hydrogen mixtures at 1 atm. 

mixtures than in the fuel-lean mixture. It is also evident from these figures that 
the ignition delay time decreases rapidly from the pure hydrogen value as the amount 
of silane is increased. Similar results were reported by Beach et al. (1980). How- 
ever, whereas the Beach mechanism predicts that the same ignition delay times will be 
obtained for temperatures above 900 K and for all percentages of added silane, the 




present mechanism predicts that over the entire range of temperature examined, the 
ignition delay time decreases as the quantity of silane increases. 

The results in figures 2(a) and 2(b) indicate that the rate of decrease of the 
ignition delay time at a given temperature becomes smaller as the mixture becomes 
richer in silane. An additional illustration of this trend appears in figures 3(a) 
and 3(b), which show ignition delay time at a given initial temperature as a function 



Silane fuel percent 



(a ) <ji = 1 . 


(b) <i> = 0.5. 


Figure 3.- Ignition delay time as a function of fuel mixture at various 
temperatures and a pressure of 1 atm. 


of silane/hydrogen fuel fraction. The increased curvature of these lines as the 
mixture approaches pure hydrogen points out the rapid sensitization of ignition pro- 
vided by adding a small amount of silane to hydrogen. Similar sensitization behavior 
has been predicted for methane ignition in the presence of ethane (Westbrook 1979). 

Examination of the time history of various reactants and products during the 
ignition period provides some further insight into the ignition characteristics of 
silane/hydrogen mixtures. A calculated time history of silane, molecular hydrogen, 
and hydroxyl radical concentrations and mixture temperature for the stoichiometric 
2 percent silane/98 percent hydrogen mixtures at constant pressure and initially at 
800 K is shown in figure 4(a). It is evident from this figure that the silane is 
consumed before any significant amounts of hydrogen begin to react. The ignition of 
the hydrogen is then enhanced (aided) by the large amounts of free radicals H, 0, 
and OH which are rapidly generated during the ignition of the silane. As the silane 
content is increased, a thermal effect (temperature increase) also contributes to the 
enhanced oxidation of the hydrogen. This effect becomes important for silane amounts 
greater than 10 percent. This combined free radical and thermal effect is evident in 




the time histories shown in figure 4(b) for a stoichiometric 20 percent silane/ 

80 percent hydrogen mixture initially at 800 K. After most of the silane has been 




Time, us 


(a) 2 percent SiH 4 /98 percent H 2 . (b) 20 percent SiH 4 /80 percent H 2 » 

Figure 4.- History of SiH 4 , H 2 , and OH concentrations and mixture temperature 
during the ignition of stoichiometric SiH 4 /H 2 mixtures at 1 atm and initially 
at 800 K. 


consumed and prior to any significant oxidation of the hydrogen, the temperature has 
increased almost 200 K. This increase in temperature and the presence of free radi- 
cals both contribute to the increased rate of hydrogen oxidation. At initial tem- 
peratures greater than 1000 K, the oxidation of silane and hydrogen is very rapid and 
takes place almost simultaneously. 

In addition to the enhancement of the ignition process through the generation of 
free radicals, silane also aids ignition to some extent by scavenging the HO 2 radical 
through the reaction sequence 


SiH 4 + H0 2 + H 2 0 2 + SiH 3 
H 2°2 + M -> OH + OH + M 


However, at temperatures above 800 K, free radical generation and the thermal effect 
are the dominant enhancement mechanisms. At temperatures below 800 K, the scavenging 
may become more important. 
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CONCLUDING REMARKS 


A chemical kinetic reaction mechanism for the combustion of silane/hydrogen has 
been assembled through analogy with methane and refined by comparing the observed 
kinetic behavior as determined in shock tube studies with that predicted by the mech- 
anism. It was determined that reasonably good agreement between the experimental and 
the calculated results could be achieved by the adjustment of the rate coefficient 
for the reaction 


SiH 3 + 0 2 ->• siH 2 0 + OH 


Calculations made with the refined reaction mechanism indicate that addition of 
small amounts of silane to hydrogen-air mixtures effectively sensitizes the mixture 
to ignition. This sensitization occurs because the silane additive rapidly oxidizes 
and provides free radicals ( H, 0, OH) in large amounts which subsequently accelerate 
the ignition of the hydrogen. As the amount of silane is increased, the rate of 
sensitization decreases. For silane amounts greater than 10 percent, a thermal 
effect (temperature increase due to the combustion of silane) also contributes to the 
sensitization effect. 

The credibility of the results given here depends to a very large extent on the 
reasonableness of the silane/hydrogen mechanism. The basic assumption of a silane 
mechanism analogous to the methane mechanism seems reasonable. However, the accuracy 
of the many assigned rate coefficients cannot be assessed. It is recommended that 
additonal kinetic studies be carried out to evaluate and improve the high temperature 
silane oxidation mechanism. More direct evidence is needed to show that the proposed 
mechanism describes the actual process. The research program should include the 
identification of intermediate species using spectroscopic techniques, matching of 
the shape of measured concentration-time profiles with those predicted, and quantita- 
tive matching of concentration-time curves with predictions. Species of particular 
interest would include SiH 3 , SiH 2 , SiH, and H. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
January 10, 1983 
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